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Multi-photon interferences with indistinguishable 
photons from independent light sources are at the 
focus of current research owing to their potential 
in optical quantum computing [H, creating re- 
mote entanglement for quantum computation and 
communication and quantum metrology 

The paradigmatic states for multi-photon inter- 
ference are the highly entangled NOON states @, 
which can be used to achieve increased resolu- 
tion in spectros cop y 05 0]? interferornetry [9i-,12l|, 
lithography 0, Il3l | and microscopy |l4l |. Multi- 
photon interferences from independent, uncor- 
related emitters can also lead to enhanced res- 
olution in metrology and imaging (li4l6l|. ho 
far, such interferences have been observ ed with 
maximally two independent emitters |l3l . 17- 26|. 
Here, we report multi-photon interferences with 
up to five independent emitters, displaying in- 
terference patterns equivalent to those of NOON 
states. Experimental results with independent 
thermal light sources confirm this NOON-like mod- 
ulation. The experiment is an extension of the 
landmark measurement by Hanbury Brown and 
Twiss who investigated intensity correlations of 
second order |l5l [. Here we go beyond this level 
by measuring spatial intensity correlations up to 
fifth order to further increase the resolution. The 
inherent simplicity of our scheme opens up the 
possibility of improved imaging in astronomy and 
biology. 

Multi-photon interference lies at the heart of optical 
quantum information processing and is an essen- 
tial ingredient in optical quantum computing P, 
quantum metrology with light and quantum imag- 
ing Multi-photon interference lies also at 
the heart of the Hanbury Brown and Twiss experiment 
[ist . One of the key signatures of classical and quantum 
multi-photon interference is super-resolution, that is, the 
appearance of fringes that are much denser than those 
obtained with standard interference techniques. 

In its simplest form the interference pattern can be 
described by the relation 



where N is the number of photons participating in the in- 
terference, Vn is the visibility, k is the difference between 
the wave vectors ki and k2 of the interfering light fields, 
with |ki| = |k2| = 27r/A, and x is the position along the 
interference pattern. When for a given numerical aper- 
ture A of the image formation system > 1, the pattern 
exhibits super-resolution and conveys information about 
source details that are smaller than the Abbe limit; when 
NA > 1 the interference pattern is sensitive to source 
structures below A [28| . This is the fundamental prin- 
ciple behind quantum lithography, where multi-photon 
interference produces sub- wavelength features [l^ ■ 

The archetypal quantum state used for demonstrating 
super-resolution is the so-called NOON state Q. Imag- 
ine a two-photon absorbing substrate that is sensitive to 
the square of the incoming intensity. If it is exposed 
to a standing light field of the form / oc cos^ ip, with 
(/? = k X, the absorption rate of the substrate scales with 



In{x) oc - [1 



VArcos(A^kx)] , 



/ oc 3 -I- 4cos(2<p) -I- cos(4(/3). If one could eliminate 
the middle term containing the slowly varying function 
cos(2(^), one would be left with an interference pattern 
displaying a fringe spacing of A/4, a factor of two im- 
provement with respect to the classical resolution limit 
Q. An iV-photon NOON state (|iV,0) + e'-^'^'\0,N))/V2 
is able to achieve this task, because the photons always 
travel together, effectively behaving like a single particle 
with a de Broglie wavelength X/N. When the two modes 
of a NOON state are made to interfere, super-resolution 
according to Eq. ^ occurs due to the A^-fold amplified 
relative phase factor cxp{iNip). However, NOON states 
are extremely difficult to realize and fragile to noise and 
decoherence; moreover, the measurement procedure re- 
quires iV-photon absorbing media. For these reasons 
NOON states are very hard to use in practice. 

Super-resolution can also be obtained with indepen- 
dent photon emitters, even though the photons no longer 
behave as a single particle with reduced de Broglie 
wavelength. In this case the slowly varying terms can 
be removed by exploiting iV"^-order intensity interfer- 
ence effects. For example, for independent photon 
sources and N detectors at carefully chosen positions 
r2,...,rAr, the iV"^-order spatial intensity correlation 
function ^'^^^(ri, . . . , r^r) as a function of ri can be 
shown to exhibit a NOON-like modulation of the form 
cos[(7V- l)(5i], with Si =^ kdsinOi (see Fig. 1) 
(1) This approach does not rely on particular states of light 
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Experimental Arrangement 



N = 2 




N = 3 




N = 4 



-•1 + 




N = 5 




Thermal Light Source gfj^g 



g^-jJsCJj.O) = j(l + -cos(Ui)| , V2!»33% 

3.5r 



-« -6 -4 -2 2 4 6 8 



ggjjis, ■"■') = ^ [l + 55 cos(2J,)] Va = 32% 
3.5|- 
3.0- 



-8 -6 -4 -2 2 4 6 



' 3 ■ 3 128 L 91 



-6 -4 -2 2 4 6 



falLSV i 2 2 3125 [ 927 J 

3.5r 




-6 -4 -2 2 4 6 



FIG. 1: iV-photon coincidence detection scheme for 
AT independent SPE or TLS (left) and calculated 
N^^ order spatial correlation functions for TLS 
(right). Left column: Detection scheme for the measme- 
ment of the normalized AT"' order spatial correlation func- 
tion g^'^\5i, . . . , 5m) for TV = 2, . . . , 5 point-like independent 
emitters located at Rq along a chain with equal spacing d 
{a = A, B, . . .) . An equal number of A'^ detectors Dj situ- 
ated at rj {j = 1, . . . , N) measure the N emitted photons in 
the far field, so that they cannot distinguish from which of 
the A' sources the photons originate. 5j = S{rj) = kdsinOj is 
the optical phase difference of two photons propagating from 
adjacent sources to detector Dj. A successful A^-photon coin- 
cident detection event is defined if all N detectors register one 
photon within a joint detection time interval. Right column: 
Theoretical plots of g^'^\5i, ...,5n) for A^ = 2, . . . , 5 TLS 
for the indicated fixed detector positions 5j {j = 2, . . . , A") 
as a function of 5i for point-like emitters (blue curve) and 
extended emitters (red curve). 



nor on A^-photon absorbing media, only on A^-fold coin- 
cidences in N separate single photon detectors. 

The occurrence of super-resolution can be used to en- 
hance the resolution in imaging and microscopy. Accord- 
ing to Abbe, an image of an object is formed if the rays 
contributing to the first diffraction order in the diffraction 
plane are captured by the imaging device since then all 
information about the object is contained in the diffrac- 
tion pattern via Fourier transform For a grating 
with A'^ slits and slit spacing d this leads to a minimal 
resolvable slit separation dmin = A/(2^), with an error 



of Adniin = X/{AA). This limit can be overcome if the 
slowly oscillating terms in the diffraction pattern of the 
grating / oc 1 + ^a=i(^ ^ ^) cos{aS) are suppressed, 
such that only the modulation at the highest frequency 
cos [(A^ — 1)6] prevails. This modulation already contains 
all relevant parameters of the grating (N and d). Based 
on counting the number of peaks M in the interference 
pattern 1 -I- Vn cos [(A^ — 1)6] across A, we find that the 
separation d and its corresponding error Ad are given by 
(see supplementary information) 



d = 



MX 



2A{N - 1) 



and 



Ad = 



A 



4:A{N-1) 



(2) 



For A^ — 1 > M > 1 the procedure allows to achieve a 
resolution beyond the Abbe limit. Here, we experimen- 
tally demonstrate an increased resolution compared to 
the canonical classical limit using up to A^ = 5. 

Consider A^ independent light sources at Rq, (a — 
A,B,.. .) along a chain with equal spacing d, and place 
A^ — 1 detectors in a semi-circle in the far field around 
the sources at specific magic angles (see Fig. 1). The 
sources have identical frequency and polarisation and 
may be single photon emitters (SPE) or thermal light 
sources (TLS). Moving another detector along the semi- 
circle and post-selecting on simultaneous single photon 
detection events in each of the A^ detectors will produce 
an interference pattern of the form of Eq. ([Ij . The inter- 
ference pattern is proportional to the (normally ordered) 
A^-point intensity correlation function 



^ {■.U,E^'\r,)E(+Hv,):) 
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(3) 



where (...) denotes the quantum mechanical expectation 
value and ^^^^^(r) and i?'^~'(r) are the positive and nega- 
tive frequency parts of the total electric field operator at 
position r, respectively. Here, £^'+)(rj) cx X]q ^''^''"^ , 
where da is the annihilation operator of a photon of 



source a and 



= |R„ — r , | is the distance between 



the source a and the detector Dj . Since the emitters are 
independent and uncorrelated, the state of the field is 
described hy p — <SS p^, where = ^„ Pa{n)\n) {n\ and 

Pa(n) denotes the photon number distribution for the 
modes originating from source a. From this, the interfer- 
ence pattern g^^") is easily calculated (see supplementary 
information) . 

The quantum paths that contribute to g^^^^ and g^^^^ 
in the case of SPE and TLS are shown in Fig. 2. For 
A^ SPE the number of quantum paths that contribute 
to the A^-photon signal is A^!, whereas for A^ TLS A^^ 
quantum paths are involved. The increased number of 
paths in case of TLS is due to the possibility of multiple 
photons originating from the same TLS, which reduces 
the visibility of g^^\ 

The form of 5*-^^ is highly dependent on the positions 
Vj of the A^ detectors. We can choose the rj such that. 
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FIG. 2: Quantum paths contributing to two- and 
three-photon joint detection events in case of N = 
2 and N = 3 emitters, respectively. Sources A and B 
are (a) two SPE [paths (I) and (II)] or (b) two TLS [paths 
(I) - (IV)]. In the case of TLS only the two indistinguishable 
quantum paths (I) and (II) contribute to the interference sig- 
nal whereas paths (III) and (IV) lead to a constant offset and 
thus a reduced visibility, (c) Sources A, B and C are three 
SPE or TLS. In case of SPE 3! = 6 quantum paths contribute 
to [first row, paths (I)-(VI)], whereas for TLS 3^ = 27 
quantum paths are involved, since in this case more than one 
photon may originate from the same source. All quantum 
paths where the same number of photons is emitted by each 
source are indistinguishable so that the corresponding quan- 
tum paths interfere coherently. The different colours refer to 
the case where one photon originates from each source (grey), 
all A'' photons originate from the same source (blue) or two 
photons originate from one source (red). The six indistin- 
guishable quantum paths in the first row in (c) are identical 
for SPE and TLS. 
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FIG. 3: Experimental setup for the multi-photon inter- 
ference measurement. The A = 2, . . . , 5 independent TLS 
are realized by opaque masks with A slits of width a — 25 fim 
and separation d = 250 fim, illuminated by pseudo-thermal 
light using a linearly polarized Nd:YAG laser at A = 532 nm 
and a rotating ground glass disc. The coherence time of the 
pseudo-thermal light source was chosen to Tc « 100 fis. The 
light from the masks is separated by 50/50 non-polarizing 
beam splitters into A beams. To measure g!^g{Si, . . . , Sn), 
A laterally displaceable fibre tips with 50 fim core diameter 
are located z ~ 1 m away from the masks, mounted on trans- 
lation stages, which collect the photons and guide them to 
single photon detectors. The output pulses of the A single 
photon detectors are fed into a coincidence detection circuit. 
GGD: ground glass disc, M: mirror, L: lens, NDF: neutral den- 
sity filter, TS: translation stage with fibre mount, BS: beam 
splitter, F: multimode fibre, Di . . . D^: photomultiplier mod- 
ules. 



and 



— cos(2(5i] 



(7) 



up to a proportionality constant, t?'-^-' takes the form of 
In-1 in Eq. ([1]), where one detector is moved and — 1 
detectors are kept at the fixed magic positions. In this 
case, for SPE, the second order correlation function in 
terms of phases Sj = kd sin 9j can be written in the form 



ffi'^(<5i,(52-0) = i[l + cos(5i)] , 
while the third order correlation function becomes 



(4) 



(3) 
9SPE 



<5i,<52 = J,<53-^) =^[l + cos(25i)] . (5) 



From Eqs. ([4]) and ([5|) one can see that for SPE the vis- 
ibility of tlie second and third order correlation function 
is 100%; this remains true for any order N . In the same 
way, the second and third order correlation function for 
TLS with the detectors at the magic positions are 



5^'L('5i,0) = - 



1 + i cos((5i) 



(6) 



displaying reduced visibilities V2 = 1/3 and V3 = 8/25, 
respectively. In all cases the super-resolving NOON-like 
modulation of g^^\5i) is clearly visible. 

Similar results are obtained for higher numbers of SPE 
and TLS. In the case of TLS, the interference pattern 
•9tls s-lways reduces to the form In-i in Eq. ([T]) when 
the detector positions are chosen as 5j = 27r(j — 2)/ {N — 
1) for J = 2, . . . ,iV. The calculated interference signals 
5tls('^i) for iV = 2, . . . , 5 independent TLS are displayed 
in Fig. 1, together with their exact analytical expressions. 

The experimental setup used to measure g^^si^^) 
shown in Fig. 3. To realize the N independent TLS, 
opaque masks with A^ identical slits are illuminated by 
pseudo-thermal light originating from a laser scattered 
by a rotating ground glass disc |30| . The large number of 
time-dependent speckles, produced by the stochastically 
interfering waves scattered from the granular surface of 
the disc, act within a given slit as many independent 
point-like sub-sources equivalent to an ordinary spatial 
incoherent thermal source. The coherence time of the 
pseudo-thermal light sources depends on the rotational 
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FIG. 4; Experimental results, (a) Measurement of aver- 
age intensities Ji and I2 at detectors Di and D2 with Di 
scanned and D2 kept constant. Since the pseudo-thermal 
light sources are spatially incoherent in first order no first- 
order interference pattern is observed, (b)-(e) Measurement 
of the normalized A^**^ order spatial correlation function 
as a function of 5i in case oi N = 2, . . . ,5 TLS for particu- 
lar detector positions S2, ■ ■ ■ , 5n ■ Red curves correspond to a 
theoretical fit taking into account the finite width of the slits. 
The only fitting parameters are the slit separation d, the slit 
width a and the visibility V}^'' . The experimentally obtained 

visibilities V^'' can be compared with the theoretical values 

(2) (3) 

Vn in Fig. 1. Single photon counting rates for grplg, g^LS' 
9tls ^^'^ 9tls about 200 kHz which with joint detection 
time windows of 50 ns, 410 ns, 410 ns and 850 ns lead to av- 
eraged Af-fold coincidence rates of 1500 Hz, 1500 Hz, 400 Hz 
and 300 Hz, respectively. 



speed of the disk and can be, for example, of the order of 
microseconds. Note that g^pLs^^^) displays the calculated 
interference signal only if the N photons are registered 
within their coherence time [l^. 

The experimental results for g^lsi^'^) 9TLsi^^) 
shown in Fig. 4. From the figure it can be seen that the 

measured curves for (7^2s('^i)' 9tls(^^) ^^'^ .9tls('^i) 
display a doubled {2di), tripled (35i) and quadrupled 
(4(5i) modulation frequency compared to 55^p^(<^i, 0) and 
f/^2^((5i,0) [see Eqs. Q and Due to the offset ap- 
pearing for TLS the visibilities Vat of g^]^g{Si) are slowly 
decreasing with N whereas for SPE the visibilities always 
remain 100%. Fig. 4 (a) exhibits the average intensities 
at detectors Di and D2 for a typical measurement of 
Stls^ demonstrating that the pseudo-thermal light used 
is indeed spatially incoherent in first order of the inten- 
sity. 



The measured curves for grls('^i) ^° 5t1s('^i) 
excellent agreement with the theoretical prediction if one 
takes into account the finite width of the slits [see red 
solid lines in Fig. 1 and Fig. 4 (b)-(e)]. The theoreti- 
cally predicted and experimentally observed visibilities. 
Vat and V^*^', are hsted in Fig. 1 and Fig. 4 (b)-(e), 
respectively. The differences between V^*^^ and Vn to- 
wards higher N are mainly due to increased dead time 
effects arising from larger joint detection time windows 
and higher single photon counting rates at the N de- 
tectors. The small deviations between the experimen- 
tal results and the theoretical curves for g^j^g and g^Ls 
are mostly due to a slight misalignment of the detector 
positions from the required magic values. From Fig. 4 
it can be seen that for a given aperture A (grey bar 
in Fig. 4) g^Lsi^'^) exhibits four times more oscillations 
than gTLsi^'^)^ which corresponds to an enhancement in 
Ad by a factor of four, and can be used to beat the Abbe 
limit. 

In conclusion, we have demonstrated spatial multi- 
photon interference patterns displaying super-resolution 
with up to five independent emitters, using pseudo- 
thermal light sources. For A'^ > 2, these experiments 
achieve a higher resolution than the canonical classical 
limit for imaging the spatial structure of the light source. 
In case of N SPE, the independent sources produce an 
interference pattern identical to the one of NOON states 
with A^ — 1 photons. The same is true for A^ TLS, except 
for a reduced visibility. Our technique does not require 
special quantum tailoring of light nor A^-photon absorb- 
ing media; it relies only on single photon detectors. The 
natural low light requirements suggest that the technique 
has potential applications for improved imaging of faint 
star clusters and in vivo biological samples. 



Supplementary information: 

First, we derive Eq. (2) in the paper. According to Abbe, 
an image of an object can only be formed if the rays 
contributing to the first diffraction order in the diffraction 
plane are captured by the imaging device . In case of 
a grating with A^ slits and slit spacing d this leads to a 
minimal resolvable slit separation and error [28 1 



TT. (8) 



"min = ^ and Admin = ^ 

where A is the numerical aperture of the imaging device 
and A the wavelength. In case of ^ = 1 we obtain dmin = 
A/2. This limit can be overcome if the slowly varying 
terms in the diffraction pattern of the grating 



/(r) cx 1 + - ^ (iV - a) cos(a<5) 



id) 



with 5 = S{r) = kd sin 9, are dispelled so that only 
the modulation at the highest oscillation frequency 
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cos [{N — 1)5] prevails. For a given A, this allows to de- 
termine a slit separation as low as dmin = X/2A{N — 1) 
with an uncertainty of Admin = X/4:A{N — 1). This can 
be seen from the following argument: let the modula- 
tion cos [{N — 1)5] display AI peaks over the range of the 
numerical aperture {—A < sinfl < A), such that 



27rA/ = 2A{N - l)kd 



MX 



2A{N -1) 



(10) 



with k = 27r/A. The uncertainty in d can be estimated 
conservatively by choosing /S.M = 1/2, what leads to 



AM 



dM 



dd 



X 



iAiN-l) 



(11) 



For N — 1 > M > 1 we can thus surpass the canonical 
classical limit given in Eq. ([8]). Counting only a single 
peak (M = 1) we obtain 



dm i n 



A 



2A{N - 1) 



(12) 



which is the new minimum resolvable slit separation. In 
case of ^ = 1 we obtain dmin = X/2{N — 1). 



J 



Next, we calculate the analytical expressions for G^^-', 
5^^) and g^'^\ For = 2 we have p = pa ® Pb 
with = a = A,B, and E^+\vj) = 

Qggifcrfjj^ With this we find 



A 



(+) 



from Eq. (3) in the paper that the intensity distribution 
G(^'(ri) at the detector Di is given by 

oo 

G'^^\yi) = ^ PA{n)PB{m) \\{n - 1, TO|i^+V, m)f 



n,m— 



+ \{n, TO — 1|_bJ'''^ m) 
{tia) + (it-b) = 2{n) ~ const , 



(13) 



where in the last equation we assumed equal amplitudes, 
that is, equal mean photon numbers (n) = {ua) = 
(aJjOa) in the modes from source a. Eq. (|13p shows that 
irrespective of the distribution Pa{n)^ no interference in 
the first order spatial intensity correlation function can 
be seen. This is due to the fact that the two sources are 
uncorrelated. 

Now, we derive g^'^\vi,Y2)- Using the result of Eq. (fT3|) 
we find 



ff^'^(ri,r2) 



1 



{{ua) + {nB)Y' 



PA{n)PB{m) l,m- 11(4+^5^+^ +B^+UJ+^)|n,TO)|2 



+ 2,TO|i*+'iJ+V,™>l'' + \{n,m ~ 2\B\'^' Bf'\n,m) 
I 



?(+)n(+)| 



(14) 



In the case where A and B correspond to single photon 
emitters (SPE) we have Pa{l) = 1 and Pa{n 1) = 0, 
so that only the first term in Eq. (fT4|) contributes to 
g^^)(ri,r2). We thus obtain 



by Bose-Einstein statistics 



P„(n) 



(l + (n„))«+i ■ 



(16) 



[1 -hcos(<52 - 5i)] 



(15) 



where (5, 



5{v,) 



k(rAj — TBj) is the relative phase 
illustrated in Fig. 1. Thus, even though the two SPE 
emit their photons independently and no information is 
contained in G^^''{ri), the normalized second order spa- 
tial intensity correlation function g^gp^(5i,62) displays a 
sinusoidal modulation with a visibility of 100% and may 
convey information about the spatial distribution of the 
sources. 

In the case where A and B correspond to thermal light 
sources (TLS), the probability distribution Pa{n) is given 



In this case, the terms in Eq. ([H]) where more than one 
photon originate from the same source do not vanish. 
The initial state |n, m) can thus lead to (^^^j^^^= 3 
different final states, namely {n — l,m — 1|, {n — 2,m| 
and (n,TO — 2|. Here, different quantum paths linking 
the same initial and final state are indistinguishable and 
therefore have to be added coherently, whereas quantum 
paths leading to different final states, being distinguish- 
able, have to be summed incoherently. From the = 4 
two-photon probability amplitudes appearing in Eq. (jl4p 
[represented by the quantum paths (I)- (IV) in Fig. 2(b)], 
only paths (I) and (II) give rise to interference whereas 
paths (III) and (IV) contribute to the offset of the two- 
photon signal. Summing the contributions for the differ- 
ent states |n, to), assuming equal mean photon numbers 
in each mode and taking into account the Bose-Einstein 
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statistics of thermal light sources Pa{n)^ we arrive at 



1 + - cos((52 - 5i) 



(17) 



Comparing Eq. (fT7|) with Eq. (fT5)) . one can see that the 
same sinusoidal modulation is obtained for two TLS as for 



two SPE, the additional terms in Eq. PT)) only reducing 
the visibility of the two-photon signal from 100% to 33%. 

For = 3 emitters, i?(+''(rj) is given by the superpo- 
sition of three fields E^+'> (r^) ^ + b" ^ 



(/('^^(ri, r2, rs) calculates to 



ikrBj 



ace 



ikrcj 



In this case 



(?(3)(ri,r2,r3) 



iinA) 



{nc)r 



P{k)P{n)P{m) 



{k - 1, n - 1, m - 1| + b(+)4+'(7(+' + . . . + Ci+'> 3^+^ A[+'>) \k, n, m) 



(fc - 2,n - l,m| 



(18) 



{k - 3, n, TO|i^+U^+Ui+^ l/c, n, m) 



j |fc, n, m) 



When ^, _B and C represent SPE, the correlation func- 
tion reduces to 

-,(3) 



.95P£;(ri,r2,r3) 
1 



27 



k{rAl+rB2+rc3) _j_ gifc(rA2+rBi+rc3) 
_|_ gifc(rAi+rB3+rc2) _j_ Qik{rA3+rB2+rci) 

_|_ gifc(rA3+I"Bl+I"C2) _|_ gife(rA2+I'S3-|-''Cl) 



(19) 



The corresponding quantum paths are shown in the first 
row of Fig. 2(c). If one places the second and third de- 
tector at positions such that 62 = 7r/4 and 63 = 77r/4, we 
obtain from Eq. (|19p 



„(3) ^ 77r 

9SPE { ^1' 4'X 



27 



[l + cos(2(5i)] 



(20) 



that is, a sinusoidal modulation which, apart from a scal- 
ing factor, oscillates at two times the frequency obtained 
for <?5p^(i5i,0) [see Eq. (|15p ]. This is identical to the 
modulation of a NOON state with N = 2. 

When the sources A, B and C represent TLS, again 
terms where two or more photons originate from the same 
source do not vanish. The initial state \k, n, m) may thus 
end up in {^^^^ = 10 different final states, leading to 

a total of N-^ = 27 different quantum paths. The corre- 
sponding three-photon transition amplitudes are shown 
in Fig. 2(c). Note that the six quantum paths listed in the 
first row in Fig. 2(c) are identical to those obtained for 
three SPE. Summing up the contributions of all possible 
quantum paths weighted with the corresponding Bose- 

J3) 



/-ON 

of three TLS to reduce dTi^si^i, 52,53) to a single sinu- 
soidal modulation with twice the frequency of ggp^((5i, 0) 
or g!^j^g{5i,0) by appropriately positioning detectors D2 
and D3. For example, if we place D2 and D3 at positions 
r2 and such that 52 ~ 0, 5^ ~ tt, we derive from Eq. 
psp . assuming again equal mean photon numbers for all 
three sources 



50 
27 



1 



25 



cos(2(5i) 



(21) 



Similar results are obtained for higher numbers of SPE or 
TLS. In the case of TLS, when positioning the detectors 
such that 

5, =2^(j-2)/(iV-l) for j=2,...,iV, (22) 

the interference pattern (7^5 always reduces to a noon- 
like modulation of the form 



9TLsi^i) '^^ + Vn cos [{N - l)5i 



(23) 



Note that the simple relation for the phase shifts 
Eq. ([221) also allows to determine A'' if the number of 
sources is unknown: according to Eq. (j22p the pure sinu- 
soidal oscillation Eq. is obtained if the TV — 1 fixed 
detectors D2 , • . . , are separated by equal amounts 
2tt/ {N~1) of their phases 52, ■ ■ ■ ,5^. Thus, by enlarging 
for different numbers of detectors N ~ 2, . . . the angles 
62, ■ ■ ■ ,0N of detectors D2, ■ ■ ■ , Dn in a uniform manner 
such that the phase relation Eq. ([221) is continuously ful- 
filled, one can monitor the interference pattern g^Lsi^^) 
until the pure sinusoidal modulation is obtained. The 
emergence of the pure sinusoidal modulation for g^Lsi^^) 
Einstein photon statistics one obtains g}p[g{5i, 52, 53). then fixes 0j, j = 2, A^, what allows to determine the 
Similar to the case of three SPE, it is also possible in case source distance d via 5j — kd sin 9j . 
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Finally, we note that the suppression of redundant data 
in order to cut out relevant spatial structural information 
as in Eq. is the goal in various schemes for super res- 
olution in microscopy, either using techniques from clas- 
sical physics (see, e.g . . IsiHsst ) or from quantum physics 
see, e.g., [H, Q, iiSi) . However, our method differs 
from these schemes as we exploit higher order {N > 2) 
non-local multi-photon interference effects, i.e., interfer- 
ence of non-local TV-photon probability amplitudes where 
the N photons appear at all N detectors at once. 
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